求解Schrdinger方程的高阶紧致差分格式

基于Richardson外推法提出了一种求解Schrdinger方程的高阶紧致差分方法.该方法首先利用二阶微商的四阶精度紧致差分逼近公式对原方程进行求解,然后利用Richardson外推技术外推一次,得到了Schrdinger方程具有O(r~4+h~4)精度的数值解.通过Fourier分析方法证明了该格式是无条件稳定的.数值实验验证了该方法的高阶精度及有效性.     

The parameter uniform numerical method for singularly perturbed parabolic reaction–diffusion problems on equidistributed grids

This article presents the study of singularly perturbed parabolic reaction–diffusion problems with boundary layers. To solve these problems, we use a modified backward Euler finite difference scheme on layer adapted nonuniform meshes at each time level. The nonuniform meshes are obtained by equidistribution of a positive monitor function, which involves the second-order spatial derivative of the singular component of the solution. The equidistributing monitor function at each time level allows us to use this technique to non-linear parabolic problems. The truncation error and the stability analysis are obtained. Parameter–uniform error estimates are derived for the numerical solution. To support the theoretical results, numerical experiments are carried out.     

A phase transition for the heights of a fragmentation tree

We provide information about the asymptotic regimes for a homogeneous fragmentation of a finite set. We establish a phase transition for the asymptotic behavior of the shattering times, defined as the first instants when all the blocks of the partition process have cardinality less than a fixed integer. Our results may be applied to the study of certain random split trees. © 2010 Wiley Periodicals, Inc. Random Struct. Alg., 39, 247‐274, 2011     

Atom waveguide and 1D optical lattice using a two-color evanescent light field around an optical micro/nano-fiber

We propose a two-color scheme of atom waveguides and one-dimensional （1D） optical lattices using evanescent wave fields of different transverse modes around an optical micro/nano-fiber. The atom guide potential can be produced when the optical fiber carries a red-detuned light with TE01 mode and a blue-detuned light with HEll mode, and the 1D optical lattice potential can be produced when the red-detuned light is transformed to the superposition of the TE01 mode and HE11 mode. The two trapping potentials can be transformed to each other for accurately controlling mode transformation for the red-detuned light. This might provide a new approach to realize flexible transition between the guiding and trapping states of atoms.     

Testifying experiment of the multi-pulse phenomena of capillary discharge soft-X-ray laser

In a capillary discharge experiment for the neon-like argon lasing, we have proposed an experimental scheme to verify that the multi-spike of X-ray diode (XRD) signal is a multi-pulse laser or is a reflection of the laser pulse in the XRD.The ceramic capillary has an inner diameter of 3mm and a length of 200mm.At the gas pressure of 28Pa and discharge current of 27kA, stable lasing has been realized.The experimental results prove that the multi-spike of XRD signal is a reflection of the electromagnetic signal produced by the laser pulse in the XRD.The improved electrocircuit scheme of the XRD to minimize the reflection phenomena is also found.     

Terahertz optical asymmetric demultiplexer based tree-net architecture for all-optical conversion scheme from binary to its other 2n radix based form

To exploit the parallelism of optics in data processing,a suitable number system and an efficient encoding/decoding scheme for handling the data are very essential.In the field of optical computing and parallel information processing,several number systems like binary,quaternary,octal,hexadecimal,etc.have been used for different arithmetic and algebraic operations.Here,we have proposed an all-optical conversion scheme from its binary to its other 2n radix based form with the help of terahertz optical asymmetric demultiplexer (TOAD) based tree-net architecture.     

The independence of initial vectors in the subdivision schemes

Starting with an initial vector λ = (λ(κ))κ∈z ∈ ep(Z), the subdivision scheme generates asequence (Snaλ)∞n=1 of vectors by the subdivision operator Saλ(κ) = ∑λ(j)a(k - 2j), k ∈ Z. j∈zSubdivision schemes play an important role in computer graphics and wavelet analysis. It is very interesting tounderstand under what conditions the sequence (Snaλ)∞n=1 converges to an Lp-function in an appropriate sense.This problem has been studied extensively. In this paper we show that the subdivision scheme converges forany initial vector in ep(Z) provided that it does for one nonzero vector in that space. Moreover, if the integertranslates of the refinable function are stable, the smoothness of the limit function corresponding to the vectorλ is also independent of λ.     

Finite Group Schemes over Bases with Low Ramification

Let A be a complete characteristic (0,p) discrete valuation ring with absolute ramification degree e and a perfect residue field. We are interested in studying the category FF _A' of finite flat commutative group schemes over A withp-power order. When e= 1, Fontaine formulated the purely linear algebra notion of a finite Honda system over A and constructed an anti-equivalence of categories betweenineFF _A'> and the category of finite Honda systems over A when p> 2. We generalize this theory to the case e – 1.     

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to non-equilibrium reacting flows is presented, Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.     

Si团簇的结构特性研究

采用由 M.Menon等人的非正交紧束缚模型,结合遗传算法和分子动力学方法,对n=2～20的中小尺寸Si团簇的结构进行了优化,计算了束缚能Eb、平均配位数CN,团簇的裂化能Efrag,总能的二阶差分Δ2E,最高占据分子轨道(HOMO)和最低未占据分子轨道(LUMO)的能隙Eg, 结合Mullikon集居数分析研究了电子特性,得出了随着团簇尺寸的增加,团簇的构型发生了近球形-长椭球形-近球形的转变、n=4,7,10等团簇具有较强稳定性、团簇中原子之间共价键性结合占支配地位等主要结论.所得的主要结论与已有的从头计算和实验结果相一致.